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We report on the observation of a unique honeycomb charge ordering of the cleaved IrTe 2 surface 
by high-resolution scanning tunneling microscopy (STM) and spectroscopy (STS). IrTe 2 was recently 
established to exhibit intriguing stripe charge orderings. Here, we show that the stripe charge order 
coexists with a metastable honeycomb phase formed locally. The atomic and electronic structures 
of the honeycomb phase are consistent with the stripe phase indicating unambiguously its charge 
order nature. A simple model of the honeycomb structure is suggested based on the overlap of three 
degenerate stripe orders, which is analogous to the 3q state description of a skyrmion. We suggest 
that the honeycomb charge order can be a route to an exotic Dirac electron system. 


The spatial ordering of charges due to the Coulomb 
interaction has been one of the key physics in com¬ 
plex electronic systems with exotic emerging quantum 
states such as high temperature superconductivity in 
cupurates [l|, [il and colossal magnetoresistance in man¬ 
ganese oxides [sl, fl- Most of long ran^ charge order¬ 
ings found so far are stripe structures [5[ while hexag¬ 
onal structures are frequently observed for two dimen¬ 
sional (2D) charge density wave (CDW) systems [1-0. 
The recent development of microscopy techniques such 
as scanning tunneling microscopy (STM) has made it 
possible to disclose atomic scale details of charge order¬ 
ings in low dimensional systems. One of the important 
achievements of the microscopic studies is the observa¬ 
tion of the checkerboard-type charge orderings in a high 
temperature superconductor 0 and a manganite per- 
ovskite liSIiil. The origin of the checkerboard order 
is still not fully understood but it corresponds to the in¬ 
trinsic or extrinsic overlap of two stripe orders |12 |. 


In this letter, we identify for the first time the hon¬ 
eycomb charge ordering for a recently found charge or¬ 
der material of IrTe2 with STM in coexistence with the 
stripe phases observed previously. The IrTe2 is currently 
under extensive investigation due to its unusual charac¬ 
teristics of the charge ordering together with the metal¬ 
lic property. IrTe2 is a layered material and exhibits 
the charge order (Ir Ir 5d^+) transition into stripe 

phases [q =(1/5, 0,-1/5)] with a large hysteresis and a 
nontrivial remnant metallic conductivity 13|-15|. This 
charge order seems to compete with the superconductiv¬ 
ity upon doping electrons, suggesting a quantum critical 
behavior 16I42Q|. The driving force of the transition is 
not fully understood yet but the CDW nature was ruled 
out [Ulii- The interplay of various different degrees of 
freedom is widely recognized, such as the strong struc¬ 
tural distortions involving the Ir-Ir dimerization and the 
depolymerization of the Te-Te interlayer bonding, the Ir 


5d charge disproportionation, the strong reorganization 
of Fermi surfaces, and the strong spin-orbit coupling. 

The honeycomb phase observed here is metastable and 
compete with the more stable stripe orders. The charge- 
ordered nature of the honeycomb structure is well sup¬ 
ported by its electronic structure revealed by scanning 
tunneling spectroscopy (STS). The microscopic model of 
the honeycomb phase is proposed based on the overlap 
of three rotationally degenerated stripe orderings, which 
is analogous to the 3q state description of a skyrmion 
and the domain wall ordering of incommensurate CDW 
systems UM- This phase has not been noticed in the 
bulk diffraction measurement and thus the surface origin 
is likely. The honeycomb charge organization in a 2D sys¬ 
tem leads readily to a Dirac electron system [25|. With 
the strong spin-orbit coupling, the honeycomb charge ar¬ 
rangement of Ir 5d states may open a route to exotic 
Dirac electron systems. 

The single crystals of IrTe2 were grown by Te flux 
using pre-sintered IrTe2 polycrystals as reported previ¬ 
ously [l3|. The sample was cleaved in a vacuum better 
than 5x10“^^ torr at room temperature or 86 K, lower 
than the transition temperature. All the STM measure¬ 
ments were obtained with a commercial ultrahigh vac¬ 
uum cryogenic STM (Specs, Germany) in the constant- 
current mode with Ptir tips at 78 K and 1.1 K. The 
differential conductance, dl/dV, was measured using the 
lock-in detection with a modulation of 1 kHz. For the to¬ 
tal energy calculation, we constructed a 6x6x1 supercell 
consisting of monolayers separated by vacuum distance 
about 2.0 nm and employed the density functional theory 
band calculation as implemented in the pseudopotential 
plane wave package, VASP 0113. For the exchange- 
correlation functional, the generalized gradient approxi¬ 
mation in PBEsol scheme has been used [28|. The spin- 
orbit effect was included. 

Figure 1(a) shows a cleaved IrTe2 surface where a Te 
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FIG. 1. (color online) STM topographies of the cleaved IrTe 2 
at 78 K. (a) A typical cleaved IrTe 2 surface shows honeycomb 
and stripe phases coexisting, (b) The FFT image of the hon¬ 
eycomb phase domain shown in (c). White circles and green 
arrows indicate the periodicities of the hexagons and atoms 
(2.5 nm and 0.393 nm, respectively), (d) The topography ob¬ 
tained at the same position as (c) after a time interval of 13 
min, showing the irreversible transition into the stripe phase. 
White arrows indicate the defect, used as the position refer¬ 
ence mark. Tunneling parameters are Vs = 1 V, T = 0.1 
nA 


layer is exposed As reported previously, vari¬ 

ous stripe phases with three different orientations (ar¬ 
rows in the figure) and with different spacings of x3, 
x5, x8, and xll are observed as inhomogeneous mix¬ 


tures p^, ll^ . We previously established that the cascade 
of stripe phases are based on the building block structures 
of x3 and x5 [Fig. 2(a)] [32|. In addition to the stripe 
phases, we find that the surface below the transition tem¬ 
perature exhibits another distinct phase with a hexagonal 
symmetry as enclosed with stripe phases of three different 
orientations. Figure 1(c) shows one such domain, which 
is composed of hexagons and honeycomb wall lattice of 
between 5 0^7 nm) and 8 times (2.7 nm) the lattice con¬ 


stant [29|, l30| . The Fourier transformation of the STM 
image [Fig. 1(b)] indicates that these hexagons have a 
quasi long range order with a rather well defined average 
periodicity of 2.5 nm. The areal ratio between the stripe 
and honeycomb phases depends on different samples, dif¬ 
ferent cleavages and on different parts of the surface but 
we always find non-negligible portion of it on the sur¬ 
face [HI. A large domain of the honeycomb phase can 



FIG. 2. (color online) Atom-resolved STM topographies of (a) 
the x5 stripe phase as well as (b) the honeycomb structure 
taken at near the Fermi energy (Vs = 20 mV). STM line 
profiles crossing the stripes and hexagons [along the dashed 
line in (b)] are also given. The white arrows indicate the 
crystalline direction of IrTe 2 . 


extend to a few hundred nm, which is readily formed by 
the rapid cooling of the sample from room temperature. 
However, the honeycomb phase is not fully stable. Very 
rarely but the honeycomb phase spontaneously converts 
to the stripe phase as shown in Figs. 1(c) and 1(d). At 
very low temperature of 1 K, the strong tunneling cur¬ 
rent pulse (in the filled state at a large bias than 0.3 V) 
can convert the honeycomb phase into the stripes or vice 
versa [HI- Otherwise, the honeycomb domains formed 
are very stable. These observations indicate that the hon¬ 
eycomb phase is a metastable phase whose energy is but 
keen to that of the stripe phases through a large energy 
barrier. 

Figure 2 shows detailed atomic scale structures of the 
stripe and the honeycomb phase. The stripe phase ex¬ 
hibits rows of atoms with the periodically modulated 
contrast, which was shown to be bias-independent and 
well characterized as due mainly to the vertical struc¬ 
tural modulation of the Te surface layer [nisi. The 
Te rows with the bright contrasts are buckled up due to 
the dimerization of Ir rows underneath [2l|, |22|. Note 
that there exist three dark rows and two bright ones in 
the most widely found stripe phase of the x5 periodic¬ 
ity, underneath which the Ir atomic rows of 5d^+ and 
5d^+ states are located, respectively. In contrast, the 
honeycomb phase shows quasi regular hexagons of the 
bright contrast as surround by honeycomb wall lattice of 
the dark contrast. However, we note that the dark and 
bright contrast are consistent in the tunneling current 
[the apparent height in the corresponding line profiles 
of Figs. 2(c)] in both phases [HI- That is, it is rather 
straightforward that the honeycomb structure consists of 
similar structural motifs, the Te buckling and ir dimers, 
to the stripe phase. 

The similarity of the stripe and the honeycomb phase 
is not limited to the atomic structure represented by the 
topographic contrast. Figure 3(a) shows the STS spec¬ 
tra taken for the bright and dark atomic rows of the 
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FIG. 3. (a) Averaged STS (dl/dV) spectra from bright and 
dark regions of both the honeycomb and the stripe phases, 
(b) A STM topography (Vs = 20 mV) of the stripe phase 
and (c), (d) and (e) spatially resolved dl/dV maps (Vs = 
-0.3 V, 0.1 V and +0.9 V, respectively) measured simulta¬ 
neously. The corresponding (f) STM topography and (g), 
(h) and (i) the dl/dV maps of the honeycomb phase. White 
dashed line indicates the single x5 stripe and a hexagon, (j) 
The height profile of the Te surface atoms (circles) of the 
monolayer hexagonal structure shown in (k), where the un¬ 
derlying Ir atoms of 4+ (red), 3+ (or less charged than 4+) 
(blue) states are also shown. The dimerized Ir atoms are con¬ 
nected by yellow rods. The 2 nm scale bar is shown in (f) and 

(j). 


stripe and the honeycomb phase. These spectra are al¬ 
most identical, indicating that the electronic structures 
of these two phases are almost the same. The prominent 
peak and dip structures in the energy range shown is 
thought to be due to the Te depolymerization or dimer¬ 
ization, which is intimately entangled with the Ir charge 
ordering and dimerization [2l|, [22|, . The charge order¬ 

ing is represented by the spectral intensity reduction in 


both filled and empty states (-0.5 - 0 and 0.3 - 1.0 V) 
on the Te rows with a bright contrast in the STM im¬ 
age. The consistent electronic structure of the stripe and 
the honeycomb phase can further be confirmed by dl/dV 
maps [Figs. 3(b) - (i)], which reveal the atomic scale lat¬ 
eral variation of the local density of states (LDOS) in 
the surface Te layer. As detailed previously [32[, in the 
stripe phase, the LDOS at +0.1 (-0.3 and +0.9) eV is 
enhanced at the bright (dark) atomic rows of the topog¬ 
raphy [^. This directly shows the charge ordering in 
the Te layer, which is coupled with that of the Ir layer 
underneath. This LDOS modulation is copied in the hon¬ 
eycomb phase as apparently shown in the figure although 
there is extra symmetry breaking in the hexagonal phase, 
the inequivalent LDOS intensity for the certain parts of 
the hexagon [Figs. 3(h) and 3(i)]. We thus can conclude 
that the stripe and the honeycomb phase have very sim¬ 
ilar atomic and electronic structure. This means that 
the bright hexagon in the STM image of the honeycomb 
phase is composed of protruded Te atoms with the dimer¬ 
ized Ir atoms underneath and the dark honeycomb 
walls of Te atoms bonded to Ir bd^^ atoms. 

The present honeycomb phase is unique among charge 
order materials while it can be compared with the 
checkerboard charge ordering for two fold symmetric 
crystals ini- In the case of CDW, domain walls of in¬ 
commensurate orderings exhibit the competition between 
the stripe and the hexagonal honeycomb phase 
Such a competition is in fact a general phenomenon 
for incommensurate superstructures, where their relative 
stability is determinded by the sign of the domain wall 
crossing energy. The uniqueness of the present system 
is that the present stripe orders are commensurate ones 
and is not a CDW system and not even insulating. 

The above similarity between the present system and 
the 2D incommensurate phases suggests that the hon¬ 
eycomb phase may be a 2D ordering in contrast to the 
strong interlayer coupling of the stripe phases In 

supporting this, the previous x-ray ion study did 

not notice this hexagonal order 13,[3l[. A different or¬ 
dering is highly likely on the surface, especially though 
the different amount of total charges on the surface layer, 
the weakening of the Te-Te interlayer bonding, and the 
possible surface strain. 

Based on the above results, we try to construct the 
structure model of the honeycomb phase from that of the 
stripe phase. As shown in Figs. 3(b), the stripe phase is 
composed of three dark contrast rows (blue atoms) and 
two bright ones (red), under which Ir^+ and Ir^+ rows ex¬ 


ist, respectively 21 M- 


Since there exist three rotation- 
ally degenerate x5 stripes, we can simply overlap these 
three structures. This idea is consistent with the fact 
the the hexagonal phase domain is always enclosed with 
stripe phases with three different orientations. This over¬ 
lap obviously yields the hexagonal unit cell with the dark 
trenches as the hexagonal walls. Within this model, one 
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can understand the varying size of the hexagons based 
on the coexistence and competition of the x3, x5 and 
x8 stripe orders. The feasibility of the present model is 
checked by the first principles total energy calculation us¬ 
ing a relatively small unit cell of 6x6x1. We found that 
the hexagonal structure is not stabilized in the bulk con¬ 
figuration but can be stabilized for the monolayer under 
the compressive strain of about 10 percent. The fully re¬ 
laxed structure is shown in Fig. 3(k), where the hexagon 
is composed of eighteen Ir dimers (yellow rods). As sug¬ 
gested above, hexagon and the honeycomb wall network 
are mainly composed of Ir'^^ and Ir^+ (or less charged 
than 4+) atoms, respectively. This hexagon corresponds 
to one of the smallest hexagons observed. The height 
distribution of the Te surface atoms in this structure is 
shown in Fig. 3(j), which is qualitatively consistent with 
the STM topography of Fig. 3(f). Further refinements of 
the model is obviously required since the LDOS maps of 
Fig. 3(g) indicates a much higher degree of honeycomb 
charge ordering than the present calculation and there is 
an extra symmetry breaking of the hexagon in the LDOS 
maps. Howeverm we suggest that this model is a reason¬ 
able starting point, which strongly indicates the charge 
order, two-dimensional and strain origin of the hexagonal 
phase. 

The essence of the present model, the microscopic 
overlap of the three rotated charge order g, is similar 
to the 3q model of the skyrmion in the 2D triangu¬ 
lar spin system [ 2 ^. The stabilization of the 3q state 
is basically a nonlinear effect due to the higher order 
term in the Landau-Ginzburg Free energy. Consistent to 
the present picture, the local-strain-induced stripe phase 
was recently reported for NbSe 2 with the host hexagonal 
CDW ordering [ 3 ^. Combined with the present finding, 
this indicates that the competition between Iq and 3q 
state is general, irrespective of the detailed mechanism 
and energetics of charge orderings. 

Very recently, the STM study of the doped samples 
with a superconduction ground state disclosed a dis¬ 
ordered honeycomb phases very similar to the present 
one [ 3 ^. We also confirmed this with our own doped 
samples. The occurrence of the honeycomb phase on the 
doped samples can be understood by the strain or the 
reduced interlayer coupling induced by dopants. This 
observation indicates the importance of the honeycomb 
charge order for the intriguing fundamental properties 
of IrTe 2 . We can further comment that the honey¬ 
comb ordering of charge in principle indicates the forma¬ 
tion of a Dirac electron system as recently demonstrated 
for the honeycomb manipulation of the surface state of 
Cu(lll) [m. The realization of such an artificial Dirac 
electronic system is under investigation currently. What 
is very much interesting in the present system is that 
such a Dirac electron system formed along the present 
honeycomb wall lattice can be very much exotic due to 
the strong spin-orbit interaction of the Ir 5d electrons 


leading to a quantum spin Hall system j38| . 

In summary, we find a metastable honeycomb broken 
symmetry phase on the surface of IrTe 2 at low tempera¬ 
ture with STM and STS. The honeycomb phase coexists 
with the stripe charge order phases identified previously. 
The structural and electronic properties of the stripe and 
the honeycomb phases are almost identical revealing the 
charge order nature of the honeycomb phase too. A pre¬ 
liminary two dimensional model of the honeycomb phase 
is suggested based on the idea of the overlap of the three 
rotationally degenerate stripe phases. The further inves¬ 
tigation of the honeycomb charge order would be highly 
interesting in creating an artificial Dirac electron system 
with a strong spin-orbit coupling. 
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